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摘 要 :基于 柱 坐 标 系 下 的 三 维 弹性 力学 基本 方程 ,采用 状态 空间 法 得 到 两 端 固 支 单 层 与 登 层 圆柱 
厚 沉 轴 对 称 问 题 的 精确 解析 解 。 为 严格 满足 固 支 端的 边界 条 件 , 将 固 支 端的 边界 位 移 函 数 作 为 状 
态 变 量 引 入 状态 方程 ,采用 增 维 方法 把 非 齐 次 状态 方程 变 为 齐 次 状态 方程 ,并 通过 层 合 渐 近 技术 将 
变 系 数 状态 矩阵 转 为 常 系数 矩阵 进行 求解 。 所 得 到 的 解 不 仅 严 格 满足 三 维 弹性 力学 基本 方程 ,而 


| 二 且 严 格 满足 国 支边 界 条 件 ,是 真正 意义 上 的 三 维 精确 解 。 算 例 表 明 ,本 研究 解 与 有 限 元 解 吻合 ,有 具 


) 有 很 高 的 精度 , 且 关于 级 数 项 数 和 分 层 数 具有 很 好 的 收 伊 性 。 另 外 ,通过 圆柱 厚 壳 各 力学 量 沿 径 向 
和 轴 向 的 精确 分 布 规律 分 析 了 厚 径 比 和 跨 径 比 变化 对 位 移 和 应 力 分 布 的 影响 
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cylindrical thick shells with two clamped ends 
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Abstract: Based on the three-dimensional basic equations of elasticity in cylindrical coordinates ,the state 
space method was used to solve axisymmetric problem of thick single and laminated cylindrical shells with 
two clamped ends. To satisfy strictly the boundary conditions of the clamped ends ,boundary displacement 
functions introduced into the state equations as state variables ,the nonhomogeneous state equations were 
transformed into the homogeneous state equations by employing the dimensional expanding method , and 
variable coefficient state matrixes were changed into constant coefficient matrixes with the approximate 
laminate technique by dividing the layer into sub-layers. The three-dimensional basic equations of elasticity 


and all the boundary conditions are satisfied strictly , the exact three-dimensional analytical solutions are 
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obtained. The results of two examples show that the present solutions obtained have high precision by com- 


paring with FEM (Finite Element Method ) solutions ,and have good constringency with the number of se- 


ries terms and sub-layers. In addition ,the influences of thickness-radius ratio and span-radius ratio on the 


distributions of the displacements and stresses of cylindrical thick shells along the radial and axial direc- 


tions are analyzed. 


Key words :clamped end;cylindrical thick shell; state equation ;boundary displacement function ;exact an- 


alytical solution 


复合 材料 圆柱 壳 因 其 几何 结构 合理 ,材料 性 能 


圆柱 壳 动力 学 问题 的 精确 解 。Santos 等 '" 采用 有 限 


优越 ,能 够 适应 于 各 种 极端 和 复杂 的 应 用 环境 ,在 航 
空 航天 、 船 舶 .机 械 和 土木 等 工程 领域 中 得 到 了 广泛 
应 用 。 复 合 材料 圆柱 壳 常 出 现 厚 度 大 、 各 向 异性 和 
不 均匀 性 等 特点 ,从 而 增加 了 力学 分 析 的 求解 难度 ; 
但 其 常 作为 结构 中 的 重要 构件 ,通常 会 提出 更 高 求 
解 精度 的 要 求 。 这 对 矛盾 使 得 传统 薄 壳 理论 不 再 满 
是 弄 解 要 求 ,因此 众多 学 者 提出 并 发 展 了 考虑 横向 
剖 玉 变形 影响 的 中 厚 壳 理论 , 主要 包括 一 阶 剪 切 理 
论 写 和 高 阶 剪 切 理论 ”) 。 为 克服 Reissner 一 阶 剪 切 
理 沦 中 前 切 因 子 中 的 影响 , Reddy 等 中 建立 了 层 合 
者 网 高 阶 剪 切 变形 理论 ,更 好 地 满足 了 过 体内 外 表 
硬 明 应 力 为 零 的 条 件 ,被 广泛 应 用 于 层 合 壳 的 求 
解 导 中 。 但 是 上 述 各 理论 均 引入 了 人 为 假设 ,三维 
弹性 力学 基本 方程 不 能 完全 相 容 , 其 求解 精度 随 着 
厚度 的 增加 开始 下 降 , 且 无 法 获得 位 移 和 应 力 沿 厚 
度 廓 向 的 精确 分 布 规律 。 
“0 为 获得 圆柱 过 真正 意义 上 的 三 维 精确 解 ,必须 
抛 挛 任何 有 关 位 移 或 应 力 模式 的 人 为 假设 ,严格 满 
是 县 维 弹性 力学 基本 方程 。 国 内 外 学 者 做 了 大 量 研 
究 久 作 '*1。 范 家 让 9 应 用 状态 空间 法 建立 了 形式 
简单 且 统 一 的 状态 方程 ,得 到 了 强 厚度 又 层 圆 柱 壳 
的 静 力 ,动力 和 稳定 性 问题 的 三 维 精确 解 , 因 状态 空 
间 法 能 很 好 地 处 理 不 同 材料 层 间 的 连续 性 问题 ,被 
广泛 应 用 于 秋 层 结构 的 精确 求解 。Tam 等 "通过 
哈密 顿 状态 空间 法 精确 分 析 了 横 观 各 向 同性 弹性 加 é 
柱 体 的 轴 对 称 问题 。 边 祖 光 等 中 采用 层 合 渐 近 模 
型 ,利用 层 间 连续 条 件 建立 了 状态 变量 的 递 推 关系 ， 
求 得 正 交 各 向 异性 功能 梯度 材料 (functionally graded 
materials, FGM) 圆柱 壳 自 由 振动 的 频率 方程 。 

传统 状态 空间 法 面 对 含 非 简 支边 界 时 通过 传 里 
叶 级 数 分 离 变 量 后 ,边界 条 件 无 法 自动 满足 ,求解 较 
为 困难 。 故 借助 数值 离散 法 来 满足 边界 条 件 的 半 解 
析 解 法 在 FGM 圆柱 厚 壳 结构 中 取得 了 广泛 的 应 用 。 
文献 [9-10] 应 用 状态 空间 法 和 微分 求 积 法 ,得 到 了 


元 法 处 理 FGM 圆柱 壳 的 轴 对 称 弯曲 问题 ,但 其 解 只 
是 近似 弹性 理论 解 。Wu 等 ' 结合 状态 空间 法 与 无 
网 格 法 ,得 到 了 含 固 支 端 FGM 夹层 空心 圆柱 的 三 维 
精确 解 。 

完全 通过 解析 手段 处 理 圆柱 厚 壳 的 非 简 支边 界 
问题 具有 特殊 的 理论 意义 和 参考 价值 ,但 现 有 文 
献 中 尚且 不 多 。 范 家 让 '“ 通过 6 函数 在 固 支 端 
引入 边界 反 力 得 到 了 固 支 丢 层 圆柱 厚 壳 的 精确 解 ， 
盛 宏 玉 等 号 则 通过 多 项 式 函 数 引 入 边界 位 移 函 数 
处 理 固 支边 界 条 件 。 刘 艳 红 等 在 此 基础 上 通过 
增 维 方法 将 非 齐 次 状态 方程 转化 为 齐 次 状态 方程 ， 
再 进行 精细 积分 ,得 到 了 相关 问题 的 三 维 精确 解 。 
但 是 上 述 解 只 是 在 分 层 点 处 满足 固 支 边界 条 件 , 并 
没有 使 固 支 端 沿 厚度 方向 严格 满足 边界 条 件 。 胡 文 
锋 等 ' 引 针对 和 矩 形 秋 层 厚 板 问 题 ,将 非 简 支边 上 的 边 
界 位 移 函 数 看 作 状 态 变量 引入 状态 方程 ,并 通过 增 
维 方法 将 非 齐 次 状态 方程 转化 为 齐 次 状态 方程 ,所 
得 到 的 解严 格 满足 固 支 边界 条 件 。 本 研究 针对 两 端 
固 支 丢 层 圆柱 厚 壳 ,为 严格 满足 固 支 边界 条 件 , 拟 采 
用 文献 [ 15 ] 的 方法 进行 分 析 , 得 到 圆柱 壳 轴 对 称 问 
题 的 精确 解 。 


1 圆柱 厚 壳 的 状态 方程 及 其 解析 解 


1.1 单 层 圆柱 厚 壳 


图 1 所 示 为 一 两 端 固 支 轴 对 称 正 交 各 向 异性 单 
层 圆 柱 厚 过 ,其 长 度 为 1, 内 径 为 RR ,外 径 为 R,, 其 厚 
度 h=R, -RR,, 中 径 ( 即 平均 半径 ) Ro = (R, + R,)/ 
2, 坐 标 轴 沿 圆柱 厚 壳 材料 弹性 主 方向 ,wv 和 w 分 
别 表示 圆柱 壳 沿 轴 向 x、 周 向 9 和 径 向 7 的 位 移 。 以 
下 从 三 维 弹 性 理论 来 分 析 过 体内 、 外 表面 受 轴 对 称 
载荷 作用 的 静 力 响应 。 


ee 


图 1 两 端 固 支 单 层 圆柱 厚 壳 
Fig.1 Single thick cylindrical shell with clamped ends 
在 柱 坐 标 系 下 , 当 不 计 体力 时 , 静 力 学 平衡 微分 
方程 为 
O00 /0% + OTg./ (7O0) +oT /or+Tr=0 
goT, wox + 90y/ (700) +oTor+2Tor =0 (1 


3 
ee 


OT,/ OX + OT%/ (7r00) + 00,/90r+(o,-0,)/r=0 
对 于 轴 对 称 问题 ,圆柱 厚 壳 的 位 移 应 力 rw、 
Tw 各 0 对 0 的 偏 导数 为 零 , 其 应 力 -位 移 关系 为 


BY [人 Co CGC 0 Ou/ Ox 

En Cs Ca Ca 0 Wi 
EI [Cs ts Ca w/or 

未 | 0 0 0 Couw[or + OW/O% 

GD 两 四 支 询 的 边界 条 件 为 

CO, x=01,u=0,w=0 (3) 


将 式 (2) 和 式 (1) 联 立 ， 消去 应 力 分 量 o, 和 ro， 
乔 训 得 到 如 下 矩阵 形式 四 , 即 


oR/or = DR (4) 
所 有 :R=[u o, 7 w]', 
0 0 Cs —Q 
| Ca/r -(Cs+l)/r -a Ca/r 
口 — Ca” Cia -1 -CQ/r ， 
Cia (C7 0 Cs/r 
a= 0/90x,C = - Ca/C3,C, = Cn - Co/Ca, C3 = 


全 人 LO CG, S03 = C0370%, CGC == C370, 
C, =1/C3 ,Cs =1/C;sso 

被 消去 的 o, 和 oo 由 式 (2) 可 以 求 出 其 形式 为 

Or， Ca -CN 0 Cr 

| [oe = 总 


为 分 离 变量 ,结合 两 端 0 现 将 各 
力学 量 展 成 如 下 形式 的 全 里 时 级 数 地 ， 


以 二 > LU(Cr)cos(CCx ) )O 二 > es ? 
= >» Tn) COS (Cx) ,Ww 二 >» W, (7) sin( Cx) 
(6) 


(5) 


其 中 :t=mm/l;m 为 非 负 整数 。 由 上 述 傅 里 叶 级 数 
形式 可 知 , 当 x=0,l 时 ,w=0 的 条 件 已 经 自动 满足 ， 
u =0 尚未 满足 , 故 将 式 (6) 中 的 w 改 为 如 下 形式 。 

Wi >» U,(r)cos(Cx) +fi(x)u'" Cr) +f(x)u? (7) 


(7) 
其 中 :w (r) 和 uw”(7) 分 别 为 固 支 端 *=0 和 x=/ 
处 的 待定 边界 位 移 函 数 ,由 边界 条 件 wu =0 确定 ; 
(x) 入 (%) 为 选 定 的 关于 x 的 函数 ,其 函数 形式 


上 可 以 设 为 多 项 式 形 式 , 故 这 里 采用 最 简单 的 
函数 , 即 
fi(x)=1-x/l, f(x)=x/l (8) 
将 函数 式 (8) 及 其 一 二 阶 导 数 也 展 成 相应 的 傅 
里 叶 级 数 , 即 


fi(x) = 1Z2 + > 4incos(Cx) 

f(x) = 172 + > Ay,cos(x) 2 

dfi(x)/dx = > Bi,sin(tx), 

dp (x)/dx = > By, sin( tx) > 

dfi(x)/dx = > Guacos(tx), 

d(x)/dx” = > Gcos( Ex) (9) 
其 中 :4，=2[1 -cos(mm)]/(mm),A,, =2[ mm *: sin 
(mm) +cos(mm) -1]/(mm),B,=2[cos(mm) - 
ll/(mml),B,, =2(1 -cosmm)/(mml),G, =0, 
G,, =06 

为 满足 边界 条 件 , 将 式 (7) 代 入 边界 条 件 u = 
令 x =0,/ ,可 以 和 号 到 
BDU) + 六 (0)z V(r) + 户 (0)u (r) = 0， 


> (U(r) +f (Du V(r) +p (Du V(r) =0 
(10) 
将 式 (6) 和 式 (7) 代 入 式 (4) 可 以 得 到 对 于 任意 
m 关于 R,, (7) 的 非 齐 次 状态 方程 , 即 
dR,(r)/dr+AW du (r)/dr+A WV du (r)/dr 
=D,R,(r) +BOu™ Cr) +BOu (r) (11) 


其 中 : 

R,(r)=[U,(r) on(r) mm(r) W,(r)], 
m_T[M2 0 0 0]7 (m=0) 

人 0 0 0]- (mz 关 0)” 


oy TI OD 0]' (m=0) 
| 0 0 0]7 (mz0)i 
0 0 Cs -¢ 
-Ci/r -(Cs+l)/r 7 Ca/r 
Ce Ci -l/r -Ce/rl 
-CC C， 0 Cs/r 
0 0 
a C3Bi,/r Be C3By,/r 
6 =0, 6, 由 - C,0, 
CIB， C1B,, 


由 于 非 齐 次 状态 方程 (11) 中 和 抢 阵 D,, .BW 和 
BO 的 系数 含有 变量 ,属于 变 系数 状态 矩阵 ,不 易 
求 枝 。 如 图 2 所 示 , 本 研究 采用 层 合 渐 近 模型 ' ,将 
音 忆 过 沿 径 向 分 成 足够 多 的 左 层 , 则 每 层 厚度 hh = 
有 对 于 任意 的 第 j (j =1,2,…,) 子 层 ,其 内 径 为 
RESR, + ( -1 加, 外 径 尺 ,= 及 + 有 中 径 Ro = 
RSP(k -0.5) 和 , 现 将 第 子 层 变 系数 状态 矩阵 
DB 和 BI 中 的 变量 > 用 第 7 层 的 中 径 Ru 表示 ， 
使 多 变 为 常 系数 矩阵 DB 和 有 0 ,形式 如 下 。 


0 0 Cs -6 
A J -Ct /Ry -(Cs+1)/R, 记 CMR 
| Cr oe = CR 


-Ct C7 0 


CR 


图 2 圆柱 厚 党 层 合 渐 近 模型 
Fig.2 Approximate sub-layered model of thick cylindrical shell 


为 求解 非 齐 次 状态 方程 (11) , 盛 宏 玉 等 "通过 
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分 层 的 方法 把 边界 位 移 函 数 wo (r) 和 zw2(r) 假 设 
成 折线 分 布 函 数 ,只 是 在 分 层 处 满足 边界 条 件 式 
(10) 。 参 考 胡 文 锋 等 "5 对 于 含 国 支边 矩形 厚 板 的 
求解 方法 ,为 严格 满足 圆柱 厚 壳 固 支 端的 边界 条 件 
式 (10) , 则 对 任意 第 7 子 层 ,该 式 关 于 r 的 导数 均 为 
零 , 即 

>» dU (7) /dr + fi(0) du (r)/dr + 


fr(0) du ® (r) /dr =0 
DC 1)"dU,(r)/dr + f(D du (r)/dr+ 


m 


(12) 


f(D du (r) /dr =0 

考虑 到 U(r) ww (r) 和 wl (7) 在 任意 两 子 层 
交界 面 均 是 连续 的 , 故 任意 子 层 均 满足 式 (12), 需 
要 圆柱 厚 壳 任 一 点 处 满足 边界 条 件 式 (10), 即 边界 
条 件 (10) 沿 厚度 严格 满足 ,这 里 以 内 表面 为 例 。 
DUR) +f (Ou (R) +fp(O)u (RR,) =0 


pd 


f(Dud (R,) =0 
(13) 
对 于 第 j 子 层 ( 即 R<r<R,,) ,将 式 (11) 和 式 
(12) 联 立 ,把 固 支 端的 边界 位 移 函 数 wo (r) 和 
wu 中 (7) 作为 状态 变量 引入 状态 方程 ,并 通过 增 维 方 
法 5 将 状态 方程 合并 成 如 下 形式 。 


OOdG(r)[dr =OG(r) (14) 
式 中 : 
Br) = [Rr) Rlr) … Ralr) wo (7) wr (r)] ， 
[TO 0 40 40 
0 TO 0 A A 
Bg FAD A 


0 0) £(0) 


0 0 

.0 0, 0, fi(l) f(D 

[DO O BY” Bi ] 

0 D, : 0 BY By 
og D B® BO | 

0 8 0 0 0 

LO 0 0 0 0 J 
其 中 :6=[1 0 0 0],0,=[(-1)"” 0 0 0], 


752 


9=[0 0 0 0],7 为 4 阶 单位 矩阵 ,O 为 4 阶 
将 状态 方程 (14) 的 两 边 同 乘 以 8; ,得 到 齐 次 
状态 方程 , 即 


dB(r) /dr =4G(7) C13 
直路 “ 9,。 
根据 状态 空间 理论 , 齐 次 状态 方程 (15 ) 的 解 为 
G@(r) =I1(r) BR,) (16) 
其 中 : 卫 (r) = ey… 称 为 状态 转移 矩阵 ;@( Rj) 为 圆 
柱 壳 第 7 子 层 内 界面 处 的 状态 向 量 。 
在 式 (16) 中 令 r= 广 ,可 得 第 / 子 层 外 界面 和 内 
界面 之 间 的 关系 式 为 
PDR, ) =1L(h) BR) (17) 
一 根据 相 邻 子 层 间 的 连续 条 件 ,从 第 1 子 层 开始 
逐 层 类 推 ,可 得 到 圆柱 厚 壳 外 表面 (第 k 子 层 的 外 表 
醒 与 内 表面 (第 1 子 层 内 表面 ) 之 间 的 关系 为 
CD G(R.) = VP(R,) (18) 
其 $= ) (hh); BR,) = 
DR,,) ;PB(R,) =D(R,)。 


1C5 又 层 圆柱 厚 壳 


四 对 于 由 疡 层 构成 的 三 层 圆柱 厚 壳 ,对 任意 的 第 i 
(人 ,2,…,p) 层 ,利用 上 述 单 层 圆柱 壳 的 方法 ,将 每 
一 拔 分 为 足够 多 的 层 ( 总 层 数 = 》 "hk; ) , 则 可 得 
到 第 i 层 的 内 界面 和 外 界面 之 间 状 态 向 量 的 关系 为 
= BR,) = BR,) (19) 
三 其 中 :R 和 R, 分 别 为 第 i 层 的 内 径 和 外 径 ; 
OD)= 了 I,(hD)…T(h4?) .1 (hM?) 为 第 i 层 的 
状态 转移 矩阵 ,hl (1<j<%) 为 第 i 层 中 第 j 子 层 的 
厚度 。 

利用 任意 相 邻 两 层 间 的 连续 条 件 

BD(R,) =B(Ri,,) (20) 

从 第 一 层 开始 逐 层 类 推 , 便 可 得 到 羞 层 圆 柱 厚 

壳 外 表面 与 内 表面 之 间 的 关系 , 即 
DR,) = YOR,) (21) 

其 中 T=. 了 , 

可 见 , 镭 层 圆柱 厚 过 的 状态 转移 方程 与 单 层 圆 
柱 厚 壳 的 形式 一 致 ,只 是 将 状态 转移 矩阵 进行 简单 
的 笃 乘 。 

通常 圆柱 厚 壳 静 力 问 题 的 内 、 外 表面 处 的 应 力 
0,(R,) 7。(R,) 和 o,(R,) rs(R.) 是 已 知 的 。 以 壳 
体内 表面 受 均 布 载荷 q 作用 为 例 ,将 内 外 表面 的 载 


答 按 式 (6) 级 数 展开 的 系数 为 
0,,(R,) =2g[L cos(mm) -11/(mm), 
Om( Ri) =Tam(R) =Tn(R) =0 (22) 

再 结合 式 (13 ) 和 式 (18) ,可 以 求 出 状态 向 量 
BP(R,) 和 BB(R,) 中 的 其 余 未 知 量 ;然后 ,从 第 一 层 
开始 逐 层 类 推 ,由 式 (17) 便 可 得 到 任意 子 层 的 
BD(R) ,再 由 式 (16) 便 可 得 到 任意 子 层 的 更 (7) ;最 
后 ,结合 式 (8) 和 式 (9) 便 可 得 到 圆柱 厚 党 任意 位 置 
相关 位 移 和 应 力 分 量 。 

将 式 (7) 和 式 (6) 中 的 第 2 式 、 第 4 式 代入 式 
(5) , 便 可 得 到 面 内 应 力 分 量 o, 和 ce 的 表达 式 , 即 
os = DE- CeU,(r) -Ciom(r) + 


CsW,(r)/rjsin(x) + Cfi'(x)u' (rr) 十 
Cfy’' (x)u'" (7) 

ou = DE- CeU,(r) -Ciom(r) + 
CaW, (r)/rl]sin(éx) + Cafi (x)u'™ (7) + 
Cf (xz)z (Cr) 


(23 ) 
由 于 固 支 端 (x =0 或 1) 处 的 9w/9r =0, 由 式 
(2) 式 (6) 式 (7) 和 式 (8) 可 得 


Or。， Cii 
中 让 u(r)] (24) 
O Ca 

2 算 例 分析 


算 例 1 一 两 端 固 支 轴 对 称 正 交 各 向 异性 单 层 圆 
柱 壳 , 内 表面 受 均 布 载荷 4 作用 , 厚 径 比 和 = PR ， 
跨 径 比 k =1R,, 取 弹性 常数 “为 

Ci,/C1, =0. 246 269 ,C1,/C,, =0. 083 1715， 

Cy/C,, =0.543 103,C;,/C,, =0. 115 017, 

C3/C,, =0.530 172,C,,/C,, =0. 159 914 

针对 厚 径 比 和 A =0.4, 跨 径 比 k=27 的 单 层 圆柱 
厚 壳 ,图 3 给 出 了 本 研究 解 与 三 维 有 限 元 (finite ele- 
ment method , FEM ) 解 的 相对 误差 随 分 层 数 上 和 级 数 
项 数 m 的 变化 关系 。 有 限 元 解 是 通过 软件 ANSYS 
采用 8 节点 的 Soild183 单元 在 精细 网 格 划分 的 情况 
下 得 到 的 结果 。 图 3(a) 表 明 各 力学 量 关 于 分 层 数 下 
的 收敛 速度 很 快 , 当 分 层 数 取 到 10 时 ,各 力学 量 
已 基本 收敛 , 固 支 端 中 径 处 切 应 力 1, 的 相对 误差 最 
大 ,为 1.7% ,其 余 位 移 和 应 力 的 相对 误差 均 小 于 


0.2% 。 图 3(b) 表 明 各 力学 量 关于 级 数 项 数 m 的 收 
伊 速度 也 比较 快 , 当 级 数 项 数 m 取 到 49 时 ,位 移 w、 
u 和 应 力 cs ae 的 相对 误差 均 小 于 0.1% ;应 力 or， 
的 相对 误差 则 波动 收敛 于 0.2% ;7 收敛 速度 相对 
较 慢 ,在 m 取 到 99 时 ,其 相对 误差 控制 在 1.7% 以 
内 。 为 了 取得 足够 精确 的 结果 ,将 取 有 =20 mm =99 
时 的 计算 结果 作为 本 研究 解 。 

表 1 给 出 了 跨 径 比 k =27, 厚 径 比 A 取 值 不 同 
时 单 层 圆 柱 厚 党 特定 位 置 位 移 和 应 力 的 本 研究 解 、 
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位 置 都 可 以 获得 十 分 精确 的 结果 ;但 本 研究 解 因 能 
够 严格 满足 固 支 边界 条 件 在 固 支 端 处 依然 能 够 给 出 
精确 的 应 力 结果 ,这 是 完全 符合 圣 维 南 原理 的 解释 ， 
固 支 边界 条 件 的 满足 情况 对 离 固 支 端 较 远 处 的 影响 
较 小 ,对 固 支 端 附近 区 域 的 影响 很 大 ;而 固 支 端 处 的 
应 力 较 大 ,往往 是 整个 结构 的 危险 点 , 故 对 其 进行 精 
确 的 计算 具有 非常 重要 的 理论 意义 和 工程 应 用 价 
值 。 需 要 注意 的 是 ,因为 有 限 元 解 在 边界 尖 角 处 很 
难得 到 精确 的 应 力 结果 , 故 表 1 中 未 取 两 端 内 外 表 


有 限 元 解 和 文献 16] 解 的 比较 。 由 表 中 可 见 ,不 同 


面 位 置 的 应 力 rc, 和 T, 结 果 。 


厚 径 比 A 的 情况 下 ,本 研究 解 和 文献 [6] 解 在 跨 中 
表 1 单 层 圆柱 过 的 位 移 和 应 力 


Tab.1 Displacements and stresses of single cylindrical shell 


， r—R, A=0.2 入 =0.4 入 =0.8 
本 本 研究 FEM 文献 [6] ”本 研究 FEM 文献 [6] ”本 研究 FEM 文献 [6] 
warh( qh) ， x=1/2 0 46.2051 46.2060 11.0267 11.0265 2.4750 2.4736 
入 0.5 44.5315 44.5295 9.9286 9.9275 1.7710 1.7700 
1 43.4489 43.4470 43.4470 9.4071 9.4060 9.4070 1.5693 1.5683 1.5690 
0 7.9584 7.9495 2.1716 2.1674 0.4757 0.4738 
0.5 2.5777 2.5782 0.6226 0.6225 0.1159 0.1158 
1 -2.7070 -2.7028 -0.8616 -0.8604 -0.1736 -0.1733 
0 2.0035 1.9675 1.9050 1.0295 1.0020 0.9660 0.5468 0.5239 0.5040 
0.5 1.7486 1.7486 0.7326 0.7325 0.2351 0.2349 
>- 1 1.5710 1.5872 1.6280 0.5705 0.5780 0.5960 0.1633 0.1656 0.1710 
vgn x=1/2 0 5.0622 5.0526 5.0390 2.6243 2.6166 2.6090 1.4827 1.4752 1.4730 
2 0.5 4.4722 4.4720 1.9542 1.9540 0.6750 0.6746 
DR 1 4.0397 4.0432 4.0530 1.5912 1.5927 1.5970 0.4551 0.4553 0.4570 
0 网 x=0 0.1 8.2398 7.9967 3.1168 3.0814 0.9694 0.9589 
_ 0.5 1.9029 1.8800 0.8351 0.8245 0.2823 0.2787 
¢ 0.9 -4.0752 -3.9674 -1.2252 -1.2131 -0.2248 -0.2221 
ob x =0 0.1 2.0292 1.9701 0.7676 0.7589 0.2387 0.2361 
0.5 0.4686 0.4629 0.2057 0.2030 0.0696 0.0686 
0.9 -1.0036 -0.9779 -0.3017 -0.2988 -0.0554 -0.0547 
Tw/g, %=0 0.1 1.6617 1.7844 1.8320 1.1431 1.1351 1.1190 0.6787 0.6681 0.6610 
0.5 1.4630 1.3831 0.7592 0.7465 0.2932 0.2916 
0.9 1.1115 1.0540 1.1390 0.5529 0.5176 0.5240 0.1815 0.1754 0.1770 


注 :本 研究 求解 时 , 取 m=99 ,k=20; 文 献 [6] 求 解 时 ,m=199,k 分 别 取 16、20 和 22。 
图 4 给 出 了 跨 径 比 x =2"7, 单 层 圆 柱 充 厚 径 比 。” 切 应 力 71 内、 外 两 面 分 布 较为 对 称 ,符合 考虑 前 切 
和 不同 取 值 时 本 研究 解 与 有 限 元 解 各 力学 量 归 一 化 ”影响 的 中 厚度 这 理论 。 当 厚 径 比 A 较 大 (和 A 宇 0.4) 
处 理 后 沿 径 向 7 的 分 布 规律 。 从 图 中 可 见 ,本 研究 ”时 ,所 有 力学 量 都 不 是 按 线性 分 布 的 , 切 应 力 7 在 
解 与 有 限 元 解 吻合 得 很 好 , 且 精 度 不 受 厚 径 比 和 A 的。 固 支 端 内 、 外 表面 附近 变化 较 快 ,位 移 w 随 径 向 变化 


影响 ,给 出 各 力学 量 的 精确 分 布 规律 ,突显 了 状态 空 
间 法 对 圆柱 厚 壳 求解 的 优越 性 。 当 厚 径 比 A 较 小 
(A<0.2) 时 , 薄 壳 的 位 移 ww 和 应 力 ws oo\ay 均 
接近 线性 分 布 , 较 符合 经 典 薄 壳 理 论 的 直 法 线 假设 ， 


较 大 ,应 力 o, (oA9 =1) 的 数值 接近 或 超过 了 应 
力 o, 和 ee 的 数值 ( 表 1) , 故 径 向 的 应 变 和 应 力 已 
不 可 忽略 ,这 说 明 传 统 的 薄 充 理论 和 中 厚度 元 理论 
不 再 适用 于 圆柱 厚 党 。 需 要 注意 的 是 ,由 于 有 限 元 
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解 在 边界 尖 角 处 的 结果 不 精确 , 故 这 里 忽略 了 1, 在 
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图 3 本 研究 解 与 有 限 元 解 的 相对 误差 随 分 层 数 和 级 数 项 数 的 变化 


Fig.3 Variation of relative error between present solution and FEM solution with the number of sub-layers and series terms 
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图 4 不 同 厚 径 比 A 下 位 移 和 应 力 沿 径 向 7 的 分 布 


Fig.4 Displacement and stress distributions with r under different thickness-radius ratios A 


图 5 给 出 了 厚 径 比 A =0.4, 跨 径 比 «不 同 取 值 
时 本 研究 解 与 有 限 元 解 部 分 力学 量 归 一 化 处 理 后 沿 
轴 向 * 的 分 布 规律 。 从 图 中 可 以 看 出 ,本 研究 解 与 
有 限 元 解 吻 合 得 很 好 ,能 够 给 出 各 力学 量 沿 轴 向 x 
的 精确 分 布 规律 。 图 5 表明 , 当 跨 径 比 k 较 小 时 ,位 
移 和 应 力 基 本 上 旦 单调 函数 分 布 ,其 端 部 效应 在 全 
轴 长 内 均 有 影响 ; 随 着 跨 径 比 « 的 增 大 , 端 部 效应 的 
影响 开始 减 小 ,只 影响 固 支 端 附近 的 区 域 , 跨 中 区 域 
趋 于 均匀 分 布 ;尤其 是 当 跨 径 比 « 较 大 (k 宇 5 ) 时 ， 
除 位 移 w 外 ,各 力学 量 逐 渐 趋 近 于 一 条 水 平 渐 近 线 ， 


该 水 平 渐进 线 正 是 6 面 内 平面 应 变 问题 ( 蜂 径 比 
K=oo 时 ) 的 解 。 由 图 5(p) 可 见 ,由 于 端 部 效应 的 影 
响 产 生 的 轴 向 位 移 v 在 端 部 有 较 大 数值 ,并 沿 轴 向 
至 跨 中 位 置 逐 渐 趋 于 零 , 故 离开 端 部 区 域 是 一 条 斜 
直线 。 本 人 研究 解 能 够 准确 地 描述 由 圣 维 南 原 理 所 忽 
略 的 端 部 效应 ,清晰 地 说 明了 路 径 比 k 对 位 移 和 应 
力 分 布 规律 的 影响 。 

算 例 2 两 端 固 支 三 层 圆 柱 厚 索 内 表面 受 均 布 载 
集 g 作 用 ,圆柱 充 长 为 1=2mRo ,厚度 为 /第 一 层 与 
第 三 层 的 厚度 hh =h =0. 1h, 第 二 层 的 厚 上 度 h, = 


0. 8j,R 为 该 三 层 壳 的 平均 半径 ,第 一 层 与 第 三 层 
的 材料 相同 ,第 二 层 的 材料 与 算 例 1 相同 ,第 一 层 与 
第 二 层 的 材料 常数 比 为 5, 三 层 厚 壳 的 各 单 层 分 层 
数 分 别 为 = =4,k, =16。 

表 2 给 出 了 路 径 比 k =27, 厚 径 比 入 取 不 同 值 
时 闭 层 圆柱 厚 壳 特定 位 置 位 移 和 应 力 的 本 研究 解 、 


有 限 元 解 和 文献 16] 解 的 比较 。E 


晶 表 2 可 见 ,本 研究 
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解 在 不 同 厚 径 比 的 情况 下 都 可 以 获得 十 分 精确 的 结 
果 , 文 献 L6] 解 仅 在 分 层 处 满足 边界 条 件 , 而 本 研究 
解 能 够 严格 满足 固 支 边界 条 件 , 故 所 得 结果 更 精确 。 
因为 有 限 元 解 在 边界 尖 角 及 交界 处 很 难得 到 精确 的 
应 力 结果 , 故 本 研究 解 两 端 内 、 外 表面 位 置 的 应 力 
cao 和 Ts 与 有 限 元 解 有 相对 较 大 的 误差 。 


表 2 又 层 圆柱 壳 的 位 移 和 应 力 


Tab.2 Displacements and stresses of laminated cylindrical shell 


ym rR, 入 =0.2 入 =0.4 入 =0.8 
本 研究 FEM ”文献 [6] ”本 研究 FEM 文献 [6] ”本 研究 FEM 文献 [6] 
wCn/( qh), =1/2 0 25.6024 25.6010 6.0408 6.0403 1.2276 1.2271 
0.1 25.4516 25.4495 5.9485 5.9478 1.1691 1.1686 
0.5 24.6023 24.6015 5.3982 5.3980 0.8736 0.8734 
0.9 23.9756 23.9735 5.0610 5.0605 0.7429 0.7427 
y lL 23.8783 23.8765 23.8770 5.0221 5.0218 5.0220 0.7326 0.7323 0.7330 
忆 0 5.8500 5.7644 5.6020 3.0114 2.9665 2.8810 1.5444 1.5128 1.4710 
= 0.1- 5.5761 5.6351 5.6490 2.7835 2.8070 2.8260 1.2900 1.3009 1.3150 
二 0.1+ 1.0307 1.0430 1.0450 0.4806 0.4856 0.4890 0.2014 0.2037 0.2070 
YY 0.5 0.9404 0.9403 0.3800 0.3800 0.1078 0.1078 
2 0.9- 0.8662 0.8590 0.8540 0.3175 0.3156 0.3130 0.0800 0.0795 0.0790 
0.9 + 4.4588 4.4253 4.3980 1.6808 1.6717 1.6570 0.4416 0.4391 0.4350 
1 4.3019 4.3432 4.4380 1.5941 1.6070 1.6440 0.4100 0.4134 0.4230 
0 14.4118 14.3900 14.3500 7.5450 7.5327 7.5140 3.9862 3.9761 3.9690 
0.1- 14.0423 14.0560 14.0600 7.1253 7.1303 7.1350 3.4174 3.4184 3.4240 
0.1+ 2.6916 2.6950 2.6950 1.3198 1.3213 1.3220 0.6053 0.6056 0.6070 
0.5 2.4331 2.4331 1.0345 1.0345 0.3197 0.3196 
0.9- 2.2297 2.2272 2.2270 0.8580 0.8580 0.8580 0.2160 0.2158 0.2160 
0.97 11.3254 11.3160 11.3100 4.4227 4.4200 4.4170 1.1379 1.1369 1.1360 
1 11.0969 11.1060 11.1300 4.2639 4.2665 4.2760 1.0688 1.0693 1.0720 
Tan/g x=0 0.1 1.8699 1.9505 1.7910 1.0462 1.2890 1.1140 0.5281 0.9308 0.6370 
0.5 1.1225 1.0746 0.5264 0.5215 0.1732 0.1725 
0.9 1.2660 1.3201 1.1290 0.4910 0.6263 0.5040 0.1214 0.2555 0.1450 
注 :”- ”和 ”+ "分 别 表示 层 间 交接 面 的 内 面 和 外 面 ;本 研究 解 mw=99, 所 + 如 +h 取 4+16+4, 文 献 [6] 解 mw=199,h +h+ 后 分别 取 3+ 


12+34+12+4 和 5+12+5。 

图 6 给 出 了 跨 径 比 x =2", 厚 径 比 A 取 值 不 同 
时 盔 层 圆柱 厚 党 本 研究 解 与 有 限 元 解 归 一 化 处 理 后 
各 力学 量 沿 径 向 7 的 分 布 规 律 。 由 图 6 可 见 , 本 研 
究 解 与 有 限 元 解 吻合 得 很 好 , 且 求解 精度 不 受 厚 径 
比 A 的 影响 ,给 出 各 力学 量 沿 径 向 的 精确 分 布 规律 ， 


近 线 性 分 布 ,位 移 w 和 应 力 o, 也 随 径 向 是 均匀 变化 
的 , 切 应 力 rs 在 内 、 外 两 端 分 布 较为 对 称 。 当 厚 径 
比 A 较 大 (A=0.4) 时 ,各 单 层 内 的 位 移 w、 应 力 or， 
和 co 不 再 按 线性 分 布 ,位 移 w 和 应 力 o, 也 随 径 向 
不 再 均匀 变化 。 应 力 o, 和 cv 在 不 同 材 料 层 间 的 分 


从 而 清楚 地 反映 出 材料 层 间 的 应 力 不 连 续 性 以 及 在 布 是 不 连续 的 ,刚度 系数 小 的 材料 中 的 应 力 要 小 于 刚 


不 同 材料 层 中 的 不 同 分 布 特点 。 当 厚 径 比 和 A 较 小 
(A<0.2) 时 ,各 单 层 内 的 位 移 v、 应 力 o, 和 re 均 接 


度 系 数 大 的 材料 中 的 应 力 , 其 他 各 力学 量 是 连续 分 
布 的 ,但 各 力学 量 在 不 同 材 料 层 分 布 规律 明显 不 同 。 
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图 5 不 同 跨 径 比 xk 下 位 移 和 应 力 沿 轴 向 x 的 分 布 


Fig.5 Displacement and stress distributions with x under different span-radius ratios K 
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图 6 不 同 厚 径 比 A 下 位 移 和 应 力 沿 径 向 7 的 分 布 


Fig.6 Displacement and stress distributions of with r under different thickness-diameter ratios A 


3 结论 


通过 将 圆柱 厚 壳 固 支 端的 边界 位 移 函 数 作 为 状 
态 变量 引入 状态 方程 ,严格 满足 了 固 支 边界 条 件 , 建 
立 友 层 圆柱 厚 壳 的 齐 次 状态 方程 ,得 到 了 释 层 圆柱 
厚 壳 轴 对 称 问题 的 精确 解析 解 。 本 研究 解 关 于 分 层 
数 和 级 数 项 数 具 有 很 好 的 收敛 性 , 当 分 层 数 上 取 10， 
级 数 项 数 m 取 49 时 ,相对 误差 基本 在 0.2% 以 内 ， 
具有 很 高 的 精度 。 该 解 能 够 给 出 单 层 与 苹 层 圆柱 厚 
壳 位 移 和 应 力 沿 径 向 的 精确 分 布 规律 , 且 其 精度 不 
受 厚 径 比 入 和 跨 径 比 « 的 影响 ,这 是 薄 壳 理论 和 中 
厚度 过 理论 所 不 能 给 出 的 。 此 外 ,该 解 因 严格 满足 
边界 条 件 , 在 固 支 端 处 依然 能 区 给 出 精确 的 应 力 分 


泊 服 现 有 精确 解法 “对 于 非 简 支边 界 的 局 
Op 
) 本 研究 方法 为 采用 状态 空间 法 分 析 一 般 边 界 条 


人 削 其 他 过 理论 或 数值 解 的 参考 依据 。 本 研究 方法 
可 帮 决 其 它 的 非 简 支 边界 条 件 问 题 ,并 不 受 材料 类 
型 条 载荷 形式 的 影响 ,具有 一 般 适 用 性 ;另外 ,本 研 
络 放 法 可 进一步 应 用 于 其 他 码 层 结构 的 热电 力 等 
多 场 看 合 的 各 种 静 ,动力 学 问题 中 。 
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